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Abstract. We report on a possible association of the recently discovered very high-energy y-ray source HESS J1825-I37 
with the pulsar wind nebula (commonly referred to as G 18.0-0.7) of the 2.1 x lO* year old Vela-like pulsar PSRBI823-13. 
HESS J1825-137 was detected with a significance of 8. 1 cr in the Galactic Plane survey conducted with the H.E.S.S. instrument 
in 2004. The centroid position of HESS J1825-137 is offset by 11' south of the pulsar position. XMM-Newton observations 
have revealed X-ray synchrotron emission of an asymmetric pulsar wind nebula extending to the south of the pulsar. We argue 
that the observed morphology and TeV spectral index suggest that HESS J1825-137 and G 18.0-0.7 may be associated: the 
lifetime of TeV emitting electrons is expected to be longer compared to the XMM-Newton X-ray emitting electrons, resulting 
in electrons from earlier epochs (when the spin-down power was larger) contributing to the present TeV flux. These electrons 
are expected to be synchrotron cooled, which explains the observed photon index of ~ 2.4, and the longer lifetime of TeV 
emitting electrons naturally explains why the TeV nebula is larger than the X-ray size. Finally, supernova remnant expansion 
into an inhomogeneous medium is expected to create reverse shocks interacting at different times with the pulsar wind nebula, 
resulting in the offset X-ray and TeV y-ray morphology. 

Key words. ISM: plerions - ISM: individual objects:PSRB1823-13, HESS J1825-137, G 18.0-0.7 - gamma-rays: observa- 
tions; 



1. Introduction 

PSRB1823-13 (also known as PSR J 1826-1 334) is a 101 ms 
evolved pulsar with a spin-down age of T = 2.1 x 10** years 
llClifton et alJll992h and in these properties very similar to 
the Vela pulsar. It is located at a distance of d = 3.9 + 0.4 
kpc (.Cor des & Lazio 2002) and ROSAT observations of this 
source with limited photon statistics revealed a compact core, 
as well as an extended diffuse nebula of size ~ 5' south-west of 
the pulsar (Finlev et al. 1998). High resolution XMM-Newton 
observations of the pulsar region confirmed this asymmetric 
shape and size of the diffuse nebula, which was hence given 
the name G 18.0-0.7 (Gaensleret al. 2003). For the compact 
core with extent /?cn ~ 30" (CN: compact nebula) immedi- 
ately surrounding the pulsar, a photon index of Fcn = 1 -6^0 2 
was measured with a luminosity of Lcn ~ '^d\ x 10^^ ergs"' 
in the 0.5 to 10 keV range for a distance of Ad^ kpc. The corre- 
sponding pulsar wind shock radius is Rs < 15" - Q.3d4 pc. The 
compact core is embedded in a region of extended diffuse emis- 
sion which is clearly one-sided, revealing a structure south of 
the pulsar, with an extension of Ren ~ 5', (EN: extended neb- 
ula) whereas the ~ 4' east-west extension is symmetric around 
the north-south axis. The spectrum of this extended component 
is softer with a photon index of Fen ~ 2.3, with a luminosity 
of Len = X lO"'^^ erg s"' for the 0.5 to 10 keV interval. No 
associated supernova remnant (SNR) has been identified yet. 

At y-ray energies, PSRB1823-13 was proposed to 
power the close-by unid entified EGRET source 3EGJ1826- 
1302 jNolan et al.ll2003h . TeV observations of this pulsar by 
the Whipple and HEGRA Collaborations resulted in only up- 
per Hmits (Hall et al. 2003; Aharonian et al., 2002.) . which are 
unconstraining with respect to H.E.S.S. 

The region around PSR B 1823-1 3 was observed as part 
of the survey of the Galactic plane with the H.E.S.S. instru- 
ment jAharonian et aljE005al) . In this survey, a source of very 
high-energy (VHE) y-rays (HESS J1825-137) 11' south of the 
pulsar was discovered with a significance of 8.1 cr. We note 
that the new VHE y-ray source is located within the 95% posi- 
tional confidence level of the EGRET source 3EG J1826-1302 
and could therefore be related to this as of yet unidentified ob- 
ject. The High Energy Stereoscopic System (H.E.S.S.) is an 
array of four imaging atmospheric Cherenkov telescopes lo- 
cated in the Khomas Highland of Namibia cHintorL2004 ). It is 
designed for the observation of astrophysical sources in the en- 
ergy range from 100 GeV to several tens of TeV. The system 
was completed in December 2003 and has already provided 
a number of significant detections of galactic y-ray sources. 
These include the first detection of spatially extended emis- 
sion fr om a pulsar wind nebul a (PWN) in very high-energy 
y-rays jAharonian et alJl2005bl). Each H.E.S S. telescope has 
a mirror area of 107 m^ (Bernlohret al. 2003) and the system 
is run in a coincidence mode ( Funketal . 2004) requiring at 
least two of the four telescopes to have triggered in each event. 
The H.E.S.S. instrument has an energy threshold of ^ 100 GeV 
at zenith, an angular resolution of ~ 0.1° per event and a point 
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source sensitivity of < 2. Ox 10 '■'cm ^s '( 1 % of the flux from 
the Crab Nebula) for a 5 cr detection in a 25 hour observation. 

2. H.E.S.S. Observations and Results 

The first H.E.S.S. observations of this region occurred as part 
of a systematic survey of the inner Galaxy from May to July 
2004 (with 4.2 hours of exposure within 2° of HESS J1825- 
137). Evidence for a VHE y-ray signal in these data triggered 
re-observations from August to September 2004 (5.1 hours). 
The mean zenith angle of the observations was 31° and the 
mean off'set (i//) of the source from the pointing direction of the 
system was 0.9°. The off'-axis sensitivity of the system derived 
from Monte-Carlo simulati ons has been confirmed via obser- 
vations of the Crab Nebula dAharonian et al.ll2005cl) . The data 
set corresponds to a total live-time of 8.4 hours after applica- 
tion of run quality selection based on weather and hardware 
conditions. 




Fig. 1. Excess map of the region close to PSRB1823-13 
(marked with a triangle) with uncorrelated bins. The best fit 
centroid of the y-ray excess is shown with error bars. The black 
dotted circle shows the best fit emission region size ((Xsoujce) as- 
suming a Gaussian brightness profile. The black contours de- 
note the X-ray emission as detected by XMM-Newton. The 95% 
confidence region (dotted white line) for the position of the 
unidentified EGRET source 3EG J 1826-1 302 is also shown. 
The system acceptance is uniform at the 20% level in a 0.6° 
radius circle around HESS J 1825-1 37. 

The standard scheme for the reconstruction o f events was 
applied to the data (see 'Aharoni an et al ]|2005d' for details). 
Cuts on the scaled width and length of images (optimized on 
y-ray simulations and off'-source data) were used to suppress 
the hadronic background. While in the standard scheme an im- 
age size cut of 80 photoelectrons (pe) was used to ensure well 
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Fig. 2. Acceptance-corrected excess slice through the region 
surrounding HESS J 1825-1 37 along the north-south direction 
of width 0.4°. The pulsar position is marked with a dotted line. 
The one-sided nature of the emission with a slow decline to the 
south of the pulsar is evident. 



reconstructed images, in the search for weak sources an addi- 
tional image size cut of 200 pe was applied to achieve optimum 
sensitivity. This cut reduces the background by a factor of 7 at 
the expense of an increased analysis threshold of 420 GeV. A 
model of the field of view acceptance, derived from off-source 
runs, is used to estimate the background. Fig.^shows an un- 
corrected excess count map of the 1 .6° x 1 .6° region around 
HESS J1825-137. A clear and extended excess is observed to 
the south of the pulsar PSRB1823-13. Assuming a radially 
symmetric Gaussian brightness profile (p oc exp(-0^ /2crl^^^^.^^)) 
an extension of (Xsource = 9.6' + 2.0' is derived. The best fit 
position for the centre of the excess lies at a distance of 11 .2' 
from PSRB1823-13 at 18h26m3s+7s, -13°45.7'+1.7'. On the 
scale of the measured displacement the systematic H.E.S.S. 
pointing uncertainties of 20" are negligible. Fig.|2lshows an ac- 
ceptance corrected excess slice through the region surrounding 
HESS J 1825-1 37, in the north-south direction, of width 0.4°. 
The position of PSR B 1823-1 3 is marked with a dotted Une. It 
can be seen that the VHE emission extends asymmetrically to 
the south of the pulsar Using the best-fit position we derive a 
statistical significance of 8.1 cr a fter accounting for all tri als in- 
volved in the search for sources lAharonian et alJ2005ah . This 
significance is obtained counting events within a circle of ra- 
dius 6 - 0.22° - 0.05 de g^), a value chosen a prio ri for the 
search for extended sources ('Aharonian et al.'"2005a'). Using a 
larger angular cut appropriate to contain the complete emission 
region of HESS J 1825-1 37 of = 0.4° an excess of 370 + 20 
and a statistical significance of 13.4 cr are derived. 

For spectral analysis, the looser image size cut of 80 pe 
and the wide angular cut of 6 < 0.4° are applied to extend the 
source spectrum to lower energies (resulting in a threshold of 
230 GeV). To reduce systematic errors, only runs with t// < 1.5° 
offset from the on-region are used, resulting in a total livetime 
of 5.4 hours. The background is estimated from regions with 
equal offset ij/ from the centre of the field of view, again to min- 
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Fig. 3. Spectral energy distribution of HESS J1825-137, as- 
suming that the X-ray emission surrounding PSR B 1823-1 3, 
the EGRET source 3EG J 1826-1 302 and the new VHE y-ray 
sour ce are related. X-ray data, indicated by lines, are taken 
from iGaensler et alJ ( l2003h and are shown for the two differ- 
ent regions as described in the text. EG RET data (full circles ) 
are taken from the third EGRET catalog ilHartman et al.ll99ql) . 
The triangles show the H.E.S.S. data from this work. 

imize systematic errors. The derived spectral energy distribu- 
tion is shown in Fig.|3] The H.E.S.S. spectrum can be fitted by 
a power law in energy with photon index 2.40 + 0.09 .s,a, + 0.2jvj 
and a flux above 230 GeV of (3.4 + 0.2„„, + 1.0,,.,) x lO"" 
cm"^s"' (corresponding to 12% of the Crab flux above that en- 
ergy). The fit has ax^/d-o.f. of 9.8/8. A spectral analysis using 
a 200 pe image cut yields consistent results. We estimate the 
systematic error on the absolute flux level to be 30%. We find 
no evidence for intra-night variability of the y-ray flux. A fit of 
the night-by-night y-ray flux of the source to a constant value 
yields a x^/d.o.f. oi 9.211. 

3. The association of HESS J1 825-1 37 with 
G 18.0-0.7 

HESS J1825-137 was discovered during an unbiased survey 
of the galactic plane within the central +30° longitude sector. 
Multiwavelength searches within a circle of radius o-^omcft ~ 
10' around the centre of gravity revealed PSRB1823-13, at 
the edge of the source radius, as the only plausible candi- 
date counterpart. The discussion below investigates the possi- 
bility of associating the pulsar and its nebula G18.0-0.7 with 
HESSJ1825-137. 

The one-sided nature of this PWN a s seen in TeV was al- 
ready suggested bv lGaensler et alJ(l20(^3h. for the X-rays, ba sed 
on the hydrodynamical simulations of'Blon din et all ( l200ll) for 
Vela X and earlier studies referenced by Gaensler et al.: it was 
assumed that the density of the medium surrounding the pro- 
genitor star was inhomogeneous along the north-south direc- 
tion, with the density towards the northern direction signifi- 
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cantly larger than to the south. The reverse shock from the 
northern direction should then have crashed relatively early 
into the PWN, pushing the latter towards the south, as ob- 
served. The apparent diameter ~ 0.5° of the TeV source in- 
dicates a relatively large PWN size, /?pwn = lld4 pc. The un- 
seen SNR shell in this scen ario would have to b e substantially 
larger: in the simulations of Blondin et alJ (l2001(). as well as in 
asample of observed composite SNRs JvanderSwaluw & Wu I 
1200 lh . the ratio /?pwn/^snr does not exceed ~ 0.25. The large 
implied Rsnr would suggest that the blast wave is expanding 
into a low-density medium to the south. For instance, a rem- 
nant in the Sedov-Taylor phase expanding into the hot phase 
of the interstellar medium, with density ~ 0.003 cm"^, would 
have^SNR = 58 pc at T = 21.5 kyears (assuming an explosion 
energy of 10^' erg). This is somewhat smaller than the value 
implied by the size of the TeV source, but uncertainties in the 
distance estimate should be kept in mind, as well as the fact that 
with a braking index different from 3 (as discussed below), the 
true pulsar age might be greater than the nominal spin-down 
time. Within the scenario outlined above, a high initial spin- 
down luminosity, or the fact that the SNR reverse shock on the 
southern side might not yet have reac hed the PWN, could also 
yield a larger value for /?pwn/^snr ( iBucciantini et alJl2004i 
and references therein). Finally, the accelerated electrons might 
diffuse beyond the boundary of the PWN, where they would 
still radiate by the inv erse Compton mechani sm but emit little 
synchrotron radiation (Ahar onian et al J 19971) . making the TeV 
source extension larger than otherwise be expected. 

The synchrotron lifetime of VHE electrons in a field of 
strength B = 10"^B_5 G, scattering cosmic microwave back- 
ground (CMBR) photons to energies Ey - lO'^fixeveV (in 
the Thomson limit) can be shown to be T(Ey) - 4.8Blg£'.pJy^ 
kyears, whereas the corresponding lifetime of keV emitting 
electrons is shorter: t(Ex) - 1 ■IB'j'^^E^^^ kyears, where /Skev 
is the synchrotron photon energy in units of keV. Gaensler et 
al. have argued that the field strength in the extended X-ray 
nebula is ~ 10 /zG, so that we assume B_5 = 1 for the X-ray 
and TeV emitting zones. The adiabatic loss and escape time 
scales from the compact nebula are a few years; this is much 
shorter than the corresponding synchrotron lifetime of X-ray 
emitting electrons, given typical values of Vela-like compact 
nebular field strengths. Synchrotron losses in the compact neb- 
ula therefore do not modify the electron spectral index of ~ 2.2 
in this region, which is responsible for the observed compact 
nebular photon index of Fcn ~ 1.6 in the 0.5 to 10 keV range. 
The shocked pulsar wind particles in the offset nebula will then 
propagate southwards until the propagation time is equal to the 
synchrotron loss time, which will result in a steepening in the 
X-ray photon index to the observed value of Fen ~ 2.3. If Vy 
and Vx are the average wind convection speeds correspond- 
ing to the respective y-ray and X-ray emitting zones, the ra- 
tio between the predicted TeV and X-ray sizes should then be 
RylRx = 4(yy/yx)B:y^(£kev/£Tev)''^ which can easily pre- 
dict a TeV nebula which is ~ 6 times larger than the X-ray neb- 
ula, given the respective experimental mean X-ray and y-ray 
photon energies of fikev ~ 2 (after absorption) and iiTeV ~ 0.9. 
This ratio (of six) also assumes that the expansion velocity of 



the PWN does not change significantly between the extended 
X-ray and TeV nebulae (i.e. Vy ~ Vx)- It should however be 
noted that this is a naive assumption and that a detailed study 
of the evolution of the velocity and associated magnetic field 
distributions should be made. 

Pulsars spin down with a braking law given by f2 = -KQ.", 
where Q. is the spin angular frequency and K is a constant, 
depending on the surface magnetic field strength and neutron 
star equation of state. Assuming typical pulsar braking indices 
n ~ 2.5 to 3, we can integrate over the spindown power E (as- 
suming a conversion efficiency of ~ 50% to electrons) to an 
epoch T(Ey) into the past when E was larger, to give the to- 
tal electron spectrum contributing to the H.E.S.S. spectral band 
(from 0.23 TeV to > 10 TeV). Since we integrate over the past 
^max ~ 98^1 kyears (corresponding to a minimum energy of 
0.23 TeV), which represents ~ 50% of the pulsar lifetime, un- 
certainties in « should still be relatively unimportant. The total 
observed TeV spectrum is then the result of a summation of 
successive inverse Compton spectra arising from past injected 
electron spectra, with the spectral break energy shifting down 
with time T into the past as T^^. The result is a cooled TeV 
spectrum, for which the TeV y-ray photon index Fjev should 
be larger than Fcn + 0.5 ~ 2.1 due to second-order Klein- 
Nishina effects on the production spectrum of IC radiation on 
the CMBR. A detailed discussion of this is beyond the scope of 
this paper, but the important point is that the observed TeV pho- 
ton index in the range 2.2 to 2.6 is consistent with this interpre- 
tation. Furthermore, because the TeV emission is the result of 
earlier epochs of pulsar injection, the ratio of y-ray luminosity 
to present spin-down power would over-estimate the true con- 
version efficiency. Compare the X-ray and y-ray energy fluxes 
in Fig. |3] the TeV energy flux is larger than the X-ray energy 
fluxes, with the latter resulting from more freshly injected elec- 
trons. 

EGRET did n ot detect pulsed emission above 100 MeV 
from this pulsar dNel et alJ Il996') and phase-resolved spec- 
troscopy is required to set upper limits on the pulsed compo- 
nent associated with PSRB1823-I3, to see if this component 
is significantly lower than the EGRET steady excess shown 
in Fig. 13 If 3EGJ1826-1302 is also associated with G 18.0- 
0.7, the GeV emitting electrons would represent the earliest 
epochs of pulsar injection. This emission should also be one- 
sided within the framework discussed above and future GLAST 
observations may be able to provide better constraints on the 
GeV morphology. 
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